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O ; Abstract 

^N I We compute the energy spectra of the gravitational signals emitted when 

a pointlike mass moves on a closed orbit around a non rotating neutron star, 

>, 

^O \ inducing a perturbation of its gravitational field and its internal structure. 

■^ 

o 
l> 

o 



c^ 



The Einstein equations and the hydrodynamical equations are perturbed and 
numerically integrated in the frequency domain. The results are compared 



Q ' with the energy spectra computed by the quadrupole formalism which assumes 

^ ' that both masses are pointlike, and accounts only for the radiation emitted 

because the orbital motion produces a time dependent quadrupole moment. 

r> ' The results of our perturbative approach show that, in general, the quadrupole 

H \ 

formalism overestimates the amount of emitted radiation, especially when the 

two masses are close. However, if the pointlike mass is allowed to move on an 

orbit so tight that the keplerian orbital frequency resonates with the frequency 

of the fundamental quasi-normal mode of the star {2ujk = ^f)j this mode can 

be excited and the emitted radiation can be considerably larger than that 

computed by the quadrupole approach. 
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I. INTRODUCTION 



The coalescence of binary systems composed of compact objects like black holes or neu- 
tron stars is considered one of the most promising sources of gravitational waves to be 
detected by ground-based interferometers. For this reason it is important to collect as much 
information as possible on the features of the gravitational signal emitted in these processes. 
This paper focuses on the phenomena which may occur during the pre-merging phase of the 
coalescence, when the two stars are still individual bodies in fast revolution around each 
other. The problem of computing the energy spectrum and the waveforms of the gravi- 
tational waves emitted in this regime can be attacked by using different approaches. The 
easiest is the quadrupole formalism, which assumes that the two stars are pointlike masses, 
and computes the emitted radiation in terms of the quadrupole moment of the system, whose 
time variation is due to the orbital motion. When the two stars get very close, post new- 
tonian corrections can be included to give a more accurate description of the trajectories, 
and to refine the orbital contribution of the emitted radiation. These calculations can be 
complemented by the inclusion of radiation reaction effects, which account for the shrinking 
of the orbit caused by the emission of gravitational waves. 

This approach clearly overlooks the fact that the two stars are extended bodies with an 
internal structure and that, when they are close, the tidal interaction becomes strong and 
new effects, unpredictable by the quadrupole -|- post-newtonian approaches, may arise. An 
accurate description of these phases of the coalescence requires the solution of the Einstein 
equations coupled with those of hydrodynamics in the non-linear regime, and many groups 
in the world are working in this direction [Q]. These studies will certainly yield new and 
interesting results, but a complete picture is still far from reaching; indeed, due to the 
complexity of this phenomenon the computational tools presently available allow to follow 
the evolution of the system for no more than a few orbits near coalescence. Waiting for 
the results of fully non linear simulations, it is interesting to explore other techniques that, 



though approximated, allow to get some insight into those phases of the coalescence where 
the quadrupole + post-newtonian approaches are inadequate. For instance, we can assume 
that one of the two stars is a "true star" , i.e. it is an extended body whose equilibrium 
structure is described by a solution of the relativistic equations of hydrostatic equilibrium, 
and that only the second star is a pointlike mass; its effect is to induce a perturbation on 
the gravitational field and on the thermodynamical structure of the extended companion, 
which can be evaluated by solving the equations of stellar perturbations in general relativity. 
In this way we can account for the tidal effects of the close interaction on one of the two 
stars and for their consequences on the gravitational emission, and get a clue on the kind 
of phenomena that could arise near coalescence. This approach has already been used in a 
previous paper p| where we computed the energy spectrum of the gravitational radiation 
emitted when a pointlike mass moves on an open orbit around a compact star. In this paper 
we shall extend our investigation to the case of closed orbits, either circular or eccentric, 
and compute the energy spectra and the waveforms of the emitted radiation. 

The purpose of our study is to compare the quadrupole radiation emitted by the sys- 
tem because of its orbital motion to the signal computed in the relativistic, perturbative 
approach. 

The orbital emission will be computed using a hybrid quadrupole approach, which as- 
sumes that the pointhke mass moves on a geodesic of the spacetime generated by the star, 
but radiates, according to the standard quadrupole formula, as if it were in flat spacetime. 

To evaluate the relativistic emission, the equations of stellar perturbations we integrate 
in the interior of the star are those derived in ref . [Q ; therefore we shall not describe them 
in detail, but only recall the relevant formulae. In p[, outside the star we integrated the 
Sasaki-Nakamura equation 0; here we consider closed orbits, i.e. a source with a compact 
support, and it is more convenient to integrate the Bardeen-Press-Teukolsky (BPT) equation 
I^J^ whose source term has a much simpler form. 

Since in this paper we do not consider the effects of radiation reaction, we cannot describe 
the evolution of the orbit and of the waveform during the inspiralling; thus the energy 



spectra we show have to be considered as representative of a certain number of orbital 
periods over which the radiation reaction effects do not produce a significant change in the 
orbital parameters of the pointlike mass (adiabatic approximation). These effects will be 
considered in a forthcoming paper. 

Stellar perturbations excited by an orbiting particle have also been considered by Ko- 
jima 0], who focused on the energy enhancement with respect to the standard newtonian 
quadrupole formula, due to the excitation of the fundamental mode by a particle in circular 
orbit. The excitation of the w-modes has been studied by Ruoff, Laguna and PuUin in a 
time domain approach [|10 . 



The plan of the paper is the following. In § II we write the equations relevant to our 
problem, in § III we discuss the source term of the BPT equation, in § IV we outline the 
integration procedure, and discuss how to find the power emitted in gravitational waves 
and the waveforms. In § V we show how to compute the same quantities by the hybrid 
quadrupole approach. Numerical results are discussed in § VI and conclusions are drawn in 
§ VII. 

II. THE PERTURBED EQUATIONS 

In order to describe the non-axisymmetric perturbations of a non rotating star induced 
by an orbiting mass we choose, as in |^ , the following gauge Q 

ds^ = e^'^'-^dt^ - e^^^^'-^rfr^ - r'dO' - r^ sin^ Odcj)^ (2.1) 

/ + 00 ( 

du e-'^'\ 2e^''NimYimdf - 2e^^' U^Yi^dr^ - 2r^H^:,de'' - 2r'' sin^ OH^idf 
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h'^^dtd<f> + h]^drdcf\ - ^^ [hLdtde + h]^drde 



^The metric (2J) in [^ contained two misprints in the terms dOdcf) and dtdcf) which have now been 
corrected. 
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The functions 



Nim, Lim, Vim, Tim 5 are the radial part of the polar (even) metric compo- 
nents, and h^^ and h}^ are the axial (odd) part. 

The unperturbed metric functions z/(r) and yU2(^) can be found by solving the 
equations of hydrostatic equilibrium (cfr. 0, Eqs. (2.2)-(2.4)) for an assigned equation of 
state. As in 0], we shall consider as a model a polytropic compact star, p = Ke"^, with 
K = 100 km'^ and n = 2. If the central density is chosen to be ec = 3 ■ 10^^ g/cm^, the 
radius and mass of the star are, respectively, Rg = 8.86 km and M = 1.266 Mq, with 
a ratio Rg/M = 4.7. Although this model is, to some extent, unrealistic, it is appropriate 
to understand the basic features of the problem we want to study. 

Inside the star, it is convenient to replace the perturbed metric functions Vim and Tim 
by two new functions 



X = nV 

77-1-1 1 U 

G = Ur[^—X -T]r + -M^' - \)\n{N + T) + AT] + ^(AT + L) - 
n r^ r 

_ e2M2 (e + p)N + -w2e2(^,-.) f^, _ T + ^^^^X] . 
2 n 

where n = {I — 1)(/ + 2)/2. The polar metric functions Nim,Lim,Xim,Gim can be 
found by solving a set of linear coupled equations (cfr. 0, Eqs. 72-75), from which the fluid 
perturbations have been eliminated 






n 



X, + —e^^'iN + L)+ cu^e^^^'-'^X = 0, 



(2.2) 



(r^G), r = nu^r{N - L) + -{e'^"^ - 1){N + L) + r(z/,, - /i2,r)A:,, + u^e^^>'^-'''^rX 



-u^^N^r = -G + Ur[X,r + ^ A^ " ^)] + " (e ^' - 1)(A^ - rX^ - r'G) 

-e^^'U + p)N + -uj^e^'^"'-''^ \N + L + -G+-[rXr + (2n + 1)X]] , 
2 [ n n ' J 

-L^r = iN + 2X).r +(-- i^,r) {-N + 3L + 2X) + 
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where Q=(|);'. 

The equations for the axial perturbation can be combined into a single wave equation by 
introducing a function Zf^{u,r) related to the axial functions by the following equations 

i d 
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(2.3) 



where r^ = Jq e "^'^'^dr. The equation Zf^{(jj,r) satisfies is (cfr. [^, Eqs. 148-149) 



drl 



Jlv 



UJ 



l{l + l)r + r^(e — p) — Qm{r) 



'ZZ = 0, 



(2.4) 



which, outside the star, automatically reduces to the Regge-Wheeler equation p. The polar 
and axial equations ( |2.2| ) and ( |2.4| ) are numerically integrated from r = 0, where we impose a 
regularity condition, up to the surface of the star r = Rg. There we compute the amplitudes 
of the Zerilli function ||^ 



zf::!{u^,Rs) 



-R.s 



nRs + 3M 



3M 



n 



Xlm{^-, Rs) — RsLlmi^, Rs) 



(2.5) 



and of the Regge-Wheeler function, Zf^{u,Rs), and their first derivatives, which will be 
used to continue the solution outside the star. 

To describe the perturbations of the Schwarzschild spacetime prevailing outside the star 
we use the BPT equation MM 



,d_ 

dr 



Adr 



+ 



{r'^uj^ + 4i{r - M)r^uj) 



iioor — 2n 



'^imi^,r) = -Tim{uj,r), (2.6) 



where A = r^ — 2Mr, and the BPT function, "^im, is related to the perturbation of the 
Weyl scalar (5\1'4 by 
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^i^{uj,r) = ^ I dQdt e^"* -2SUe,<f)) \r^ 6^^{t,r,e,<f)) 
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where 



2Slm{(^, ' 



is the spin-weighted spherical harmonic 
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The advantage of using the BPT equation is that (5\l/4 is invariant under gauge transfor- 
mations and infinitesimal tetrad rotations, and that the squares of its real and imaginary 
parts are proportional to the energy flux of the outgoing radiation at infinity in the two 
polarizations. In addition, the source term T/m(co',r), which will be discussed in detail in 
the next section, has a much simpler form than the source of the Sasaki-Nakamura equation 
which we used to study the gravitational emission in the case of open orbits in . 
At the surface of the star, where Tim{uj, r) = 0, the relation between "^im, 



Zf^ and 



^Im 
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Im 
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According to ( p.9|) we can write \E'/m as a sum of an axial and a polar part, i.e. 

vl/, = vl/,'^^ + vl> P"' 

^ im ^ Im ^ ^ Im ^ 
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However, in Section 4 we will show that depending on the value of the harmonic indices I 
and m, only the polar or the axial part of "^im have to be selected. 

III. THE SOURCE TERM OF THE BPT EQUATION 

We shall assume that the pointlike mass mo which excites the perturbations of the star 
follows a geodesic of the unperturbed spacetime on the equatorial plane, with energy E 
and angular momentum L^, so that the geodesic equations are 

dt E dr n^ / 2M\ ( Z^V dip L^ , , 

For a closed orbit the motion takes place between a periastron rp and an apoastron r^, 

roots of the equation 7 = (a third root marks the onset of a plunging motion, that is not 

relevant for our study). We can define the semi-latus rectum p and the eccentricity e in 

terms of rp and ta through the relations: 

pM pM 

^^ = TTi' '^ = T^e ^'-'^ 

Both "p and e are dimensionless, and they are, respectively, a measure of the size and of the 

degree of noncircularity of the orbit. Note that < e < 1. The orbit is periodic in the radial 

coordinate, and quasi-periodic in the yj-coordinate, i.e. 

r(t + At)=r(t) (3.3) 

if>{t + At) = ip>{t) + Av5. 

The source term of the BPT equation ( |2.6| ) is 

I T^ 

Tim{uj, r) = -2\Jn{n + l)r''T(„)(„)/„(u;, r) - 2^/n^K+—T(^n){m)imiyJ, r) (3.4) 

A r^ 
- ^A+— A+rT(^)(^)/m(u;,r) , 



and can be found as follows. The stress-energy tensor of the orbiting mass 

T- = E ^ ^ ^ ^(* - t'.ir)y'' i» - f ^('■)) . (3.5) 

where tk{r), ^k{r) are the time and angular position of the mass on the kth semi-orbit, 
is projected onto the Newman-Penrose tetrad, {l,n,m,m), to find its tetrad components 
T{p){q) = T'^'^e(p)^e(q)^. These are subsequently expanded in the spin- weighted spherical 
harmonics sSim{(^,(p), and Fourier expanded, as follows 

T(p)(,)/™ {u, r) = ^J dtdQe'-' ^S^ i^) T^p^ (t, r, Q), (3.6) 

where s = — 2, — 1, for T(^rn)(m),T(n)(m),T(n)(n) respectively. As a result of this procedure 
we find 

\k=—oo j 

The explicit expressions of Ti^(^p')i^q)im are given in appendix A. By making use of the relation 

+ 00 +CXD 

Y^ e*^^ = 277 Y. ^ (^ - 27rj) , (3.8) 

A:=— oo j=—oo 

it is easy to see that Eq. (p.7|) can be written as 

27r +°° 

T{p){q)lm {Uj,r) = — Y ^{^- ^mj) To(^p)iq)lm (^, r) , (3.9) 

j = -co 



where 



^rnj = -r- = J^lr + mil^, (3.10) 

and the source of the BPT equation takes the form 

+00 

Tira{uJ,r)= Y S{UJ -UJmj)To lm{,UJ,r). (3.11) 

j=-ao 

It should be noted that Qr and Q^ are the two characteristic frequencies of the problem. The 
frequency fi^ = -^ is associated to the periodicity of the radial motion, whereas ^ip = -^ 
is the angular velocity of an inertial observer with respect to which the yj-motion of the 
orbiting mass appears to be periodic. 
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IV. THE SOLUTION OF THE BPT EQUATION 



The inhomogeneous BPT equation ( |2.6| ) can be integrated by constructing a Green func- 
tion which ensures that (5\l/4 matches regularly with the interior solution at the boundary 
of the star, and behaves as a pure outgoing wave at infinity. This problem has been solved 



by Detweiler in the case of black holes [|rT|; here we shortly describe the simple general- 
ization of the method to the case of stars. First we discuss some symmetry properties 
of the functions involved in the problem we want to solve. Under complex conjugation, 
spherical harmonics behave as Y';J^(6', 0) = (— l)'"yi_m(6', 0); consequently, the perturbed 
metric functions in ( p.l|) and the functions Zf^ and Zf^, satisfy the following property 



Fi*^{r,uj) = {-l)"'Fi_m{r, -oo). From eqs. (piP it follows that 



i^Ci^,r)r = i-ir^' ^i-^'mi-^^r) (4.1) 

By inspection of the source term, we find that \E';m must satisfy an additional relation 

v^,:,(r,^) = (-l)'vl>,_^(r,-.;). (4.2) 



In order ( |4.1| ) and ( [4.2| ) be consistent, and looking at eqs. ( |2.1CI| ) we see that the following 
selection rule must hold 



if (/ + m) is even, "^im = "^ 



pol 
Im ' 



• if (/ + m) is odd, ^i^ = ^i^\ 

Thus, depending on the value of the harmonic indices / and m, "^im is either polar or axial. 
As explained in Section II, we integrate the equations of stellar perturbations in the interior 
of the star ( p.2|) and (p^), and construct the functions Zj^^{uj,Rs) and Zf^{uj,Rs) and 
their first derivatives at r = Rg] from them we compute '^i^{uj,Rs) and '^i^ (u;,i?s), 
as given in eqs. ( p.ll[ ), and their first derivatives, which are needed to integrate the BPT 



equation outside the star. However, it should be noted that the regularity condition imposed 
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at r = allows to determine Z^^ and Zf^ only up to an unknown amplitude, Ximi^), 
to be determined by the matching conditions at the boundary of the star. In what follows, 
we shall indicate as ^i^{uj, Rs) and ^i^ {u, Rg) the values of the axial and polar part 
of the wavefunction "^im as computed by numerical integration of the interior equations. 
The problem we want to solve therefore is 



Ij rpt"^] mi uj,r) = -T,m.(uj,r 



lm\ 



(4.3) 



*j„(r ^ oo) = r^e'^''*Aim(uj) 



"^'imi^^Rs) = Xlm{uj)^'lmi^,Rs), 

where IL^pr is the differential operator on the left hand side of the BPT equation. If [l+m) 
is even, ^ im{uj , R^) = ^ i^\uj , Rs) , whereas if (l + m) is odd ^ im{uj , R^) = '^ i^ {u , Rs) . 
Aim{uj) is the unknown wave amplitude to be determined. The general solution of eqs. ( |4.3|) 



IS 



'^im(uJ,r) 



w. 



Im 



^/m. / -TTT ^/m. Tim + '^Irr, I "T^ '^ Im Tim 



A 2 ^ Im -'- Im ~r ^ im 

Rs '-i 



(4.4) 



where \E';^ and \E';^ are two independent solutions of the homogeneous BPT equation 
defined as 



lLBPT^,^(^,r) = 



'^ilni^,r~^oo)=r'e 



S'lujr, 



flLsPT^,^(^,r) = 

'^ll{uJ,Rs)=^lm{uJ,Rs) 
['^lmi^,Rs) = ^l'm{^,Rs 



(4.5) 



and Wim{i-o) is the Wronskian 
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yyimK^J . ^Im^lm ,r ^ Im^ 



Im ,r 



(4.6) 



From eq. 



it is easy to see that the amplitude of the wave at infinity is 

dr' 



Alm.iuj) 



WirrdUj) JRs A^ 



(4.7) 



By the use of eq. (|3.11| ) this expression becomes 
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1 l-K ^°° t°° dr' 

AiU^) = - ... , . -TT E -^ (^ - ^-i) / Ti ^«m(^' '^O ^0 ;m(^, r'). (4.8) 

H/jml^j At ^.^_^ JKs A^ 

Some further details related to the evaluation of the amplitude ( [4. 81 ) are given in Appendix 
A. 

We shall now compute the time-averaged energy-flux 

/^\=lto^=lin,ii:/dJ^) , (4.9) 



where the energy spectrum, ^, can be expressed in terms of the wave amplitude at infinity 



Aim{uj) as 

Since the wave amplitude can be written as (cfr. Eq. |4.8| ) 

oo 
Aim{uj)= Y^ Aimiuj)S{uJ - UJmj) (4.11) 



j = -oo 



we have 



1 +°° r 1 

^li^y E E Jduj"6{u-u^,r ^,\Aim{^)\' (4.12) 



'rf^\ ,. 1^+2?/- ,„ 1 

— ) = iim — 

, dt / T^oo T 



-.i? 



Im j=—oo ' 

+ 00 

= E E , \ |Am(^mj)P = E E ^Lj, 
im j=-oo ^'"^mj im i=-oo 

where "S"^" is the regularized squared (5-function, such that 

27r 



Iim -^"6{u — uJmjY" = S{^ — ^mj), (4-13) 



and we have defined the time-averaged power spectrum 



ELj = -r^\Am{^mj)r. (4.14) 

In conclusion, the gravitational emission is characterized by a series of spectral lines at 
frequencies Umj- From the symmetry properties 

^-j-m = ~^mj (4-15) 
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and 

A*^{umj) = i-iy ii-„(a;_,_„), J > (4.16) 

it follows that 

Ei_^_j = Ei^j. (4.17) 

Thus, once we know the power spectrum Ej^j as a function of the frequencies Umj, for an 
assigned value of / and for positive m, the spectrum for negative m is obtained by Eq. ( [4.171) . 



Since (5\E'4(t, r, ^, 0) and the gravitational wave amplitude in the radiation gauge are 
related by 

5<il,{t, r, e, </>) = -^ [hl'^it, r, e, 0) + i/^7(t, r, 9, 0)] (4.18) 

using Eq. ( p.7|) we find 

rhl'',^{t,r,9,4>)] = _2St^{9,0) Ree'^^'f' ^ ^ i,„(,;„^.)e— -(*"-) (4.19) 

J r— >oo . (^ 

j=— oo rraj 

oo r) 

r/i7,^ (t, r, 0, 0)1 = _25,„ (0, 0) Im e*'"'^ ^ ^ i,„(^^,)e-*--(*-^*) 

j=-oo "^mj 

where we have separated the 0— dependence of the spin weighted spherical harmonics, 

-2Slm{0A)=-2Slm{9,0)e'^t 

V. THE QUADRUPOLE EMISSION 

We shall now compute the gravitational radiation emitted by the mass mo because of its 
accelerated orbital motion around the star. The energy flux is computed by using a semi- 
relativistic approximation, which assumes that mo moves along a geodesic of the curved 
spacetime, but radiates as if it were in flat spacetime. Using the quadrupole formula, it is 
easy to show that the TT-components of the gravitational wave emitted by the particle are 
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rh'J{t,r,e,<j)) 



+ 



[Qxx - Qyy) COS^ + Qxj, sill 20 
^yy ^ zz cos (7 -|- L^ zz '^xx 



(1+ 



cos'^ 



(5.1) 



rh^^{t, r, 6^, 0) = - cos 6 (Q^^ - Qyy^ sin 20 + 2(5 ^^ (^ ~ 2 cos^ 



where Qki denotes the components of the reduced quadrupole moment, Qm = 
rrio (X^X' — i^f |X| j, and 0,0 are the polar angles. The two-dimensional vec- 
tor X is the position of the particle along its trajectory in the equatorial plane X = 
{r(t) cos Lp(t),r{t) simp (t)), and r(t) and (p(t) are given by the geodesic equations 
( p. 11) . The expressions of the second time-derivative of the components of Q^i in terms of 
r{t) and ip{t) are 



Qxx = fno (a cos 2ip — j3 sin 2ip + 5/2) , 

Qyy = niQ (—a cos 2ip + P sin 2ip + 6/2) , 

Qxy = mo {a sin 2ip + jS cos 2ip) , Qzz = ^moS, 



(5.2) 



where 



a = r'^ + rr — 2r^(y9^, 



(5.3) 



P = Arr^p + r^(/3, 5 = 2 (r^ + rf j /3. 



In Appendix B we explicitely compute the Fourier transform of the metric components ( p.lj) , 
which will be used to evaluate the energy flux, and we show that they can be written as 



h^J{uJ, r,9,(f)) = Yl Yl ^(^ ~ ^rnj) Hee^Umj, 0, 

m=— 2,0,2 j=~oo 



(5.4) 



V (^'^'^' 



m=—2,2 j=—oo 



where uj^j are defined in eq. ( |3.10D , and Hgg,Hg^ are given in Appendix B. We shall 
now derive the time-averaged quadrupole energy flux (|4.9|). Since 
rfE(Q) 1 



dSdt IGtt 



hig'{t,r,e,<f)) + hiut,r,e,(t>) 



(5.5) 



it follows that 
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using Parseval's theorem this becomes 

,2 oo 



/dE^Q^\ _ r 
\ dt / ~ 16^ 



where 
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m=-2,2' 



E E ^: 



mj5 



j=-oo m=-2,0,2 



for m = —2, 2 



for 771 = 



"^^ 167r 



t/fi 



c/fi 






^mj \He 



(5.7) 



(5.8) 



VI. NUMERICAL RESULTS 



The equations of stellar perturbations ( |2.2| ), ( |2.4| ) and ( p.6| ) have been numerically in- 
tegrated for a set of bounded orbits identified by selected values of the orbital parameters 
{E, Lz), or, equivalently, (e, rp). In computing the energy flux, we have seen that the energy 
is emitted at a discrete, infinite set of frequencies Umj, with — cxo < j < +oo, defined in 
eq. ( p.lOl) . The output of our perturbative calculations are the amplitudes of the spectral 
lines El^j (EiH) and the corresponding waveforms (|4.19|) . 

The energy computed by the hybrid quadrupole approach is also emitted at the same 
discrete frequencies Umj', however, whereas the quadrupole emission is resticted to Z = 2, 
and m = (—2,0,2) for hj^ and m = (—2,2) for hJJ (cfr. Eqs. ^.4|) , for the relativistic 
calculations / > 2 and —l<m<l for both polarizations. Thus to compare the outcome of 
the two approaches we have to confront the quadrupole spectral lines E^, with the / = 2 



-mj 



relativistic lines Ef^-. 

In Figure |I] we show, as an example, the energy output for an orbit with periastron 
rp = 3Rs and eccentricity e = 0.1 computed by the quadrupole approach, and the relativistic 
results for Z = 2, 3 and 4, for the same orbit. The spectral lines are plotted for the discrete 
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values of the dimensionless frequency Mumj, for assigned values of positive m. We do 
not plot the lines corresponding to negative m because they can be obtained through a 
reflection across the zero frequency axis of the positive ones, by virtue of the symmetry 
property ( [4.17 ). A comparison of the quadrupole emission (upper panel, left) with the / = 2 



relativistic emission (upper panel, right), shows that for 171 = and m = 2 the two spectra 
are qualitatively similar. As expected, the three plots which refer to the perturbative results 
show that most of the energy is emitted in the / = 2 multipole, and that for each /, the 
/ = m component is always larger than the others. 

It is known that for particles in circular orbit around black holes the total power emitted 
in each multipole 

oo 2 

El'= i: E EL, (6.1) 

j=—oo m=—2 

scales with the multipole order as 

Ef^ ~ p^-', (6.2) 



where p is the orbital semi-latus rectum |T^ . This power-law scaling was found analytically. 
Subsequently, Cutler, Kennefick and Poisson [|l5l numerically integrated the BPT equation 
for a Schwarzschild black hole with a point particle moving on bounded orbits. They showed 
that the same result holds, at least in order of magnitude, also for particles in eccentric orbits. 
We find that a similar power law exists also for stars, both for circular and eccentric orbits, 
as indicated in Figure where, as an example, we plot the ratio E^ / E2 as a function of 
/, for an orbit with e = 0.1 and p = 15.64 {rp = 3 Rs). In Table 1 we tabulate the ratio 
E^ / E2 for different values of / and for two circular orbits with tq = 3 i?^ and tq = 10 i?^, 
respectively. 

In Figure ^ we show how the energy output obtained by the relativistic approach varies 
as a function of the eccentricity of the orbit; we consider four cases, e = 10~^, 10~^, 0.1, 0.4. 
All orbits have the same periastron {rp = 3Rs), and the plots are given for / = rra = 2, since 
this is the dominant contribution to the emitted radiation. In the zero eccentricity limit, 
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the whole power is concentrated in the harmonic with j = 0, corresponding to a frequency 
Ucirc = 2-^ = 2c<jfc, where Uk is the keplerian orbital frequency; as the eccentricity increases, 
the frequency of the highest line slightly decreases, and higher order harmonics become 
significant. 

We have compared the total power computed by the quadrupole formalism 

oo 

i^^= E E <• (6-3) 

j=-oo m=-2,0,2 

with that emitted in the / = 2 multipole, E2, defined in eq. (|6.1|) and computed by 



the perturbative approach, for different values of the eccentricity and of the periastron. We 
find that, in general, E2 is sistematically smaller than E'^. The amount of emitted 
radiation affects the orbital evolution of the system and the shape of the gravitational signal, 
in particular during the latest phases of coalescence, where the orbit is already circularized. 
To understand the relevance of this effect, which may be important for the detection of 
these signals by the ground based interferometers VIRGO and LIGO, we have computed 
the relative difference 

-^-. 6.4 

EQ ^ ' 

for circular orbits, as a function of the orbital radius, tq. The results are shown in Figure ^. 
For large values of the radius the relative difference tends, as expected, to zero; at a distance 
of 10 stellar radii it is about 7 %, and it becomes greater than 14 % when the two stars are 
3 Rs apart. 

In order to check the correctness of our results, we have repeated the calculation by a 
different approach, integrating the inhomogeneous Zerilli and the Regge- Wheeler equations 
for the same orbits. The results agree to the round-off error. 

The situation changes if the point mass moves on an orbit "resonant" with a mode of 
the star, which means the following. For the model of star we are considering, the lowest 
frequency mode is the fundamental one, whose frequency is UfM = 0.12034. To excite 
this mode the mass should move on an orbit such that the frequency of one of the spectral 
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lines of the quadrupole emission, ujmj, with m = —2, 0, 2, is very close to, or coincides 
with Uf. We find that, as the quadrupole spectral line frequency uj^j approaches 
ujf for some value of m and j, the amplitude of the emitted radiation, computed in the 
perturbative approach, increases. This suggests that the excitation mechanism could be seen 
as a resonant scattering of the gravitational wave emitted by the system in the orbital motion 
(the quadrupole wave) on the potential barrier generated by the perturbed star. Indeed, 
the discrete nature of the power spectrum emitted in the quasi periodic motion of the point 
mass suggests an analogy with an atomic laser: in this picture, the atomic energy levels 
correspond to the quasi-normal mode frequencies, and the quadrupole radiation frequencies 
to the energy of the electromagnetic radiation exciting them. 

In order to excite the fundamental mode of our star, the two bodies must be very close, 
and therefore it is reasonable to assume that the orbit is circular and that Umj reduces to 
ujcirc defined above. To show how efficient this resonant mechanism could be, in Figure ^ we 
plot the energy output E22J, and the corresponding quadrupole energy E2J for a point 
mass moving on an orbit such that uj^rc = ^f {^o = 1.37417 Rg). From this figure we see 
that the situation changes dramatically with respect to the non-resonant case: the energy 
emitted in the relativistic calculation is about 600 times larger than that computed by the 
quadrupole approximation. Whether the fundamental mode could be excited during the 
coalescence of neutron star binary systems is, however, questionable, and will be discussed 
in the concluding remarks. 

It should be mentioned that the /-mode excitation by a particle in circular orbit around 
a star was studied also by Kojima 0. In his paper he only considered circular orbits and 
polar perturbations with 1 = 2, m = ±2. We find the same qualitative behaviour, although 
with minor differences of the order of 10% in the wave amplitude. We are confident on 
the correctness of our results since, as mentioned above, we got them using two completely 
different formalisms, one based on the Regge- Wheeler, the other on the Newman-Penrose 
approach. 

In Figure P, we show the + polarization of the waveform - the x polarization is zero 



because we assume the observer is on the equatorial plane - for a circular orbit with radius 
tq = 3Rs, and for / = 2. The gravitational waveform obtained in the quadrupole ap- 
proximation is also shown (dashed lines) for comparison. There are two remarkable effects. 
The first is that the m = 1 axial contribution to the relativistic waveform, which is usually 
ignored, is not negligible. Indeed it induces a beating of the axial and polar frequencies, 
clearly seen in the figure. The second effect is that the average of the amplitudes of the 
positive (or of the negative) peaks in the relativistic waveform, is smaller than that of the 
quadrupole waveform. This is related to the fact that the / = 2 relativistic energy output 
is systematically smaller than that of the quadrupole as we get close to the star (see Figure 
3). 

A case with large eccentricity (e = 0.4) is shown in Figure 0. The structure of the 
waveforms is now much more complicated. However, both effects seen in the circular case are 
still present. The beating of the axial and polar frequencies produces similar changes in the 
maxima of the wave amplitude, i.e about 10 % in both cases. We have found that the relative 
contribution between the axial and polar emission El^ ""^ /^i ^° is quite independent of 
eccentricity and decreases approximately as ~ 1/p. Thus, it becomes negligible at large 
distances but it might be significant in the late stages of the inspiralling. 

VII. CONCLUDING REMARKS 

In this paper we have studied the gravitational emission of a binary system composed of 
a star and a point mass orbiting around it, by using a perturbative approach. The results 
have been compared with the orbital emission computed by the quadrupole formalism, which 
assumes that both objects are pointlike, thus neglecting the fact that the star is an extended 
body with an internal structure, and that the dynamical evolution of the gravitational field 
couples to the thermodynamical evolution of the star. Of course the perturbative approach 
also is quite a crude approximation of a realistic binary system, since one of the two stars 
is still considered as a point mass. However, it allows to treat at least one star in an 



19 



exact manner, since its internal structure and its gravitational field are exact solutions of 
the equations of hydrostatic equilibrium. The interaction with the companion is treated 
as a perturbation, and evaluated by linearizing the Einstein equations coupled with the 
hydrodynamical equations. 

In comparing the outcome of the two approaches, we find the following effects: 

1) The total relativistic power E2 emitted in the / = 2 multipole is smaller than that 
computed by the quadrupole approach, E'^ . 

2) The waveforms have a different shape, due to the fact that the m = 1 axial contribution, 
which is absent in the quadrupole scheme, produces a beating of the axial and polar fre- 
quencies. 

3) If the point mass is allowed to get sufficiently close as to excite the fundamental mode 
of the star, the amplitude of the wave computed by the relativistic approach significantly 
increases with respect to the quadrupole prediction. About this point, it should be noted 
that the frequency of the fundamental mode of neutron stars, Uf, is expected to be of 
the order of ~ 2 — 3 kHz, depending on the equation of state prevailing in the interior. 
This frequency is too high to be excited, since the coalescing system would reach the ISCO 
(Innermost Stable Circular Orbit) and merge, before the resonant orbit is reached. However, 
ujf is affected by rotation, and for fast rotating stars it may become small enough to be 
"excitable" . In that case the emitted radiation would probably be enhanced by the rotation, 
and this is an interesting effect that has never been studied in a relativistic framework. 

In this paper we have considered a very simple model of star described by a polytropic 
equation of state. More realistic models of neutron stars allow the existence of other classes 
of modes at lower frequency with respect to Uf. For instance, the (7-modes may lay in 
a frequency range of about ~ 100 — 400 Hz, which is the region where the ground-based 
interferometers are more sensitive. If the resonant excitation of these modes is efficient, 
the signal emitted during coalescence may suddenly change when the resonant orbit is ap- 
proached, introducing a new feature in the expected waveform. 

Besides this, several other issues remain to be clarified. The waveforms we produce by the 
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relativistic approach are different from those computed by the hybrid quadrupole approach. 
Although this is more accurate than the newtonian quadrupole approach, because it assumes 
that the pointlike mass moves on a geodesic of the unperturbed spacetime, we are aware of 
the fact that by the PPN formalism (see e.g. [0) it is possible to refine the trajectory of the 
particle especially near coalescence, and have a more accurate evaluation of the radiation 
emitted because of the orbital motion. Thus the question is: is the difference between the 
signal we compute and the most accurate estimate of the signal which is provided by the 
PPN formalism still significant? We believe the answer is positive, because the difference 
between the relativistic and the quadrupole signals we find can be attributed to the role 
played by the internal structure of the star, and to the way in which the gravitational field 
couples with the fluid. However, this question has to be answered by a direct comparison. 

In order to produce waveforms that can be used as templates in the data analysis of 
gravitational wave experiments, radiation reaction effects have to be considered. We are 



working to include these effects in our scheme following refs. p5| , p!6| , and to evaluate how 
the evolution of the system changes with respect to the traditional picture. Important 
questions that need to be answered to construct a matched filter and to extract the chirp 
mass of the coalescing system are: i) what is the number of cycles the gravitational signal 
does in the bandwidth of the interferometers, ii) how the amplitude changes in time, iii) 
how much these effects depend on the equation of state of dense matter. All these issues 
will be considered in subsequent papers. 
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APPENDIX A: THE SOURCE TERM OF THE BPT EQUATION 

In this appendix we discuss the procedure to find the wave amphtude Aim{uj) (cfr. 

Eq-El) 

1 27r +^ _ . r"^ dr' , 

lit -L 



Am{^) 



Wl,r,(uj) 



— J2 H^- ^mj) J^ -^ *«m(^, r) To imiuj, r') 



(Al) 



j=-oo 



where 



/ T /\ T 

To im{r, uj) = -2^n{n + l)r%(„)(„) - 2^/n^^+—To{^n){rh) - —K+^^+rTo(rh){rh) ■ (A2) 
The tetrad components of the stress energy tensor of the pointhke mass are 



To(n){n)lmi!^jf) 



I ( E^\ ^ , . 2E . , 
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^L, 



7r^ 



— sin(cc;t — mys) H cos(cc;t — my?) 

cos(a;t - mip) _2S't^{'K/2, 0) , 



0'S'L(vr/2,0) 
i5;j7r/2,0), 



where r, t(r) and ^p{r) refer to the point mass trajectory. From the expression of Tq im{f^i ^) 
we see that the last two terms contain the differential operator A+ applied to a function 
of r; when we evaluate the integral (|A1|) these terms can be integrated by parts defining the 



(A3) 



operators K± = ^A±, where A± = -^ ±iuj, and using the property 



Rs 



dr f {r)k+g{r) 



Rs 



drg{r)K_f{r) 



which holds if, as always in our case, /(r) or g{r) vanishes at the extrema of the integration. 
After applying this procedure, and replacing the expressions of the To(p)(q)im^ the wave 
amplitude can be written as 

1 27r +°° 



Almi'^) 



Wi„,iu) 



— Yl ^{^ - ^mj) J dr [Jo + /_i + /_2] 



(A4) 



j=~oo 



where 
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Since the integrand of Eq. (^1) diverges at the turning points where 7 = 0, it is convenient 
to perform the numerical integration using a different integration variable, defined by the 
following equation 

■pM 



"<-{x) 



(A5) 



1 + e cos X 
where x ranges from to 27r in a whole orbit. In terms of x the equations of motion become: 



dt 



p'^M{p - 2 - 2e)i/2(p - 2 + 2e)i 



/2 



dx (p ~ 2 — 2ecosx)(l + ecosx)2(j9 — 6 — 2ecosx)^''^ 
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p 



1/2 



d(p 
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APPENDIX B: FOURIER TRANSFORM OF THE QUADRUPOLE WAVE 



In this appendix we explicitely compute the Fourier transform of the metric components 
of the wave emitted by the pointlike particle because of its orbital motion, given in Eq. 
( p. ID , and discuss the symmetry properties of the corresponding spectrum. The orbit is 
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quasi-periodic in (y9, but exactly periodic in (r, r, f, 0, 0) so that we can decompose the 
Fourier transform as a sum over periods 



h{uj) 



1 



(Bl) 



(■(n+l/2)At 
„=-oo A«-l/2)Ai 

where (to(^), V^ol'")) indicates the branch of the trajectory which starts at the periastron 
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h{tn = h + nAt, (pn = ^0 + nA(p)e "dtn, 



(r = rp) at t = and ends at the apoastron (r = r^) at t = ^; (— to(^), — V^o(^)) indicates 
the "mirror" branch starting at ta and ending at rp, with t G [— -g^, 0]. An inspection of Eqs. 



( ^.1|) and ( p.2|) shows that the integrals are essentially of three types: 1) those containing 
6, say hi(t), which do not depend on ip and therefore are exactly periodic; 2) those that 
contain a periodic term, say h2{t), times sin2(y9, and 3) those that contain a periodic term, 
say hslt), times cos2(/?. These integrals can be developed in the following way: 
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To obtain this result we have used Eq. ( p.8| ) and the definition ( p.lO| ) of uj^j. Similarly we 
have: 
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By using this procedure we find that the wave components can be written as 



where 



and 
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In terms of these integrals, the various contributions to the average power radiated 
appearing in formula (|5.7| ) are given by 

,2 



■/3 

2iy ' 



i5;2^ = ^(^2,)^(X2",+X^ 



■/3 



2iy ' 



0, ,2 /r'5^ 



Using Eq. (|4.15| ) and the definitions of the various integrals, it is straightforward to prove 
that: 



X' 



X' 



2-j -^2i' 



xt 



2-i 



^2i5 



2-i - 

^2-i 



J' 



-2i' 



2i' 



/C^ 



/C^ 



(B3) 



25 



Finally, from Eqs. ( [4.15| ) and (P3|) we find 



E'l^_,=E% (m = 2,0,-2) (B4) 



that are similar to the properties ( 4.17 ) for the relativistic energy fiux 



Thus the quadrupole power spectrum has essentially the same frequency content and 
the same symmetry properties as the the relativistic spectrum with / = 2 and m = —2, 0, 2, 
except for the m = 1 contribution, which is missing. 
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TABLES 
TABLE L Ratio of the power in the ^-th multipole, E[^, to the power in / = 2 multipole for 
two circular orbits of radius tq = SRg and tq = WRs, respectively. 



El'/E, 



R 



3Rq lORq 



3 8.6-10-2 2.7-10-2 

4 9.2 - 10"^ 9.2 - 10-"^ 

5 1.0 - 10-3 3.3 - 10-5 

6 1.2 -10-"^ 1.2-10-6 
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FIG. 1. The gravitational emission associated to an eccentric orbit ■with e = 0.1 and rp = 3Rs, 
is illustrated by plotting the amplitude of the spectral lines versus the dimensionless frequency 
Mujmj- The radiation computed by the hybrid quadrupole approach E^- (upper panel, left) 
has to be compared 'with that computed by the relativistic approach, EI^-, for I = 2 (upper 
panel, right). In the lower panel 'we plot the relativistic I = 3 and I = 4 contributions. 
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FIG. 2. The ratio of the total power emitted in the /-th multipole, E[^, to the total power 
emitted in / = 2, £"2 , is plotted as a function of I, for an orbit with e = 0.1 and p = 15.64 
{rp = 3 Rs). It clearly shows a power law behaviour E[^ ~ p'~^. 
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FIG. 3. In this figure we show how the spectral content of the gravitational emission changes 
as a function of the eccentricity, plotting the spectral lines Elj:^- as a function of the dimensionless 
frequency Mujmj, for I = m = 2 and for the same value of the periastron rp = 3 Rs- As the 
eccentricity increases, the location of the highest line shifts slightly towards a lower frequency, and 
higher order harmonics become more relevant. 
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FIG. 4. The relative difference between the total power E^ computed by the hybrid 
quadrupole approximation, and the total power emitted in the I = 2 multipole, E2 , computed 
by the relativistic approach, is plotted for circular orbits as a function of the radius tq (given in 
units of the stellar radius). When the point mass moves on an orbit far from the star the two 
approaches give the same result; the quadrupole emission becomes significantly larger than £"2 
for ro <WRs. 
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FIG. 5. The spectral line emitted by the system when the point mass moves on a close 
circular orbit is compared to the same quantity computed by the hybrid quadrupole formalism. 
The relativistic line (in black) is much larger than the quadrupole one (in white), because the 
frequency of the quadrupole line coincides with that of the fundamental mode of the star, ujf, 
and a mechanism of resonant excitation occurs (see text). 
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FIG. 6. The /i^ component of the gravitational wave emitted when the point mass moves on 
a circular orbit with vq = 3 Rs, is plotted versus the retarded time in units of the orbital period. 
Since we assume that the observer is on the equatorial plane, the h^ component vanishes. In 
order to compare the relativistic waveform (continuous line) with the waveform computed by the 
hybrid quadrupole approach (dashed line), only the / = 2 component of the relativistic signal is 
shown. The difference between the two signals is basically due to the m = 1 contribution of the 
axial perturbations to the relativistic waveform (see text). 
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FIG. 7. As in figure 6, we plot the /iq^ component of the gravitational wave emitted when 
the point mass moves on an eccentric orbit with e = 0.4 and rp = 3Rs- The structure of the 
waveform is now much more complicated, but the beating of the frequencies induced by the m = 1 
axial contribution is still present. 
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